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The heating of carriers in an intrinsic graphene under an abrupt switching off a dc electric field 
is examined taking into account both the energy relaxation via acoustic and optic phonons and the 
interband generation-recombination processes. The later are caused by the interband transitions 
due to optical phonon modes and thermal radiation. Description of the temporal and steady-state 
responses, including the nonequilibrium concentration and energy as well as the current-voltage 
characteristics, is performed. At room temperature, a nearly-linear current-voltage characteristic 
and a slowly-varied concentration take place for fields up to ~20 kV/cm. Since a predominant 
recombination of high-energy carriers due to optical phonon emission at low temperatures, a de- 
pletion of concentration takes place below ~250 K. For lower temperatures the current tends to be 
saturated and a negative differential conductivity appears below ~170 K in the region of fields ~10 
V/cm. 

PACS numbers: 72.80.Vp, 72.20.Ht 



I. INTRODUCTION 

In addition to detailed investigation of the linear trans- 
port phenomena in graphene, see reviews [l| and refer- 
ences therein, a high-held transport regime is coming un- 
der examination during recent years. Starting from the 
first experiments on current cleaning of graphene, Q a set 
of measurements on transistor structures, where lateral 
inhomogeneities of concentration and contact phenom- 
ena are essential, were performed, see last results and 
reviews. [3, H{ Hot electron transport in short-channel 
suspended graphene devices was studied in Ref. 5. Re- 
cently, a few measurements [H, 0] were performed on 
homogeneous structures without contacts contribution. 
As for theoretical description, the energy relaxation pro- 
cesses in graphene were studied in several papers. [8] 
Heating of carriers by a strong dc electric field was con- 
sidered both analytically d, |Io[ and numerically, with 
the use of the Monte Carlo method. [ll| Most of these 
results and HUSE! were obtained for the monopo- 
lar graphene, with a fixed concentration of electrons (or 
holes). For the case of intrinsic graphene, not only the 
energy (temperature) of carriers increases due to the 
Joule heating but also a carrier concentration changes. 
Thus, the nonequilibrium distribution is determined 
both the energy relaxation and the interband generation- 
recombination processes. This regime of heating was an- 
alyzed in [9j for the low-energy carriers, at electric fields 
E <10 V/cm and at low temperatures. The modifica- 
tions of the current-voltage characteristics, I(V) (here V 
is a drop of voltage along structure), with variation of 
the gate voltage, V g , from the heavily-doped regime of 
transport to the intrinsic case (at V g =0), were reported 
in Q for the room temperature. While the mechanism 
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for formation of the second Ohmic law in I(V) under 
the monopolar regime of transport was considered in [7| , 
an appearance of the quasilinear characteristic at V g =0, 
when the generation-recombination processes are essen- 
tial, remains unclear. 

In this paper, we consider a heating of carriers after an 
abrupt switching off a dc electric field. In addition to the 
quasielastic relaxation via acoustic phonons and the in- 
terband transitions due to thermal radiation considered 
in Ref. 9, the intra- and interband emission and absorp- 
tion of optical phonons is taken into account in the high- 
energy region. The Cauchy problem for the quasiclassical 
kinetic equation is solved below for the case of a weakly 
anisotropic distribution under an effective momentum re- 
laxation. The temporal concentration and energy as well 
as the current-voltage characteristics are analyzed and 
their dependencies on temperature and field strength are 
described for the steady-state conditions. 

The results obtained can be briefly summarized as fol- 
lows. At room temperature, a nearly-linear I(V) char- 
acteristic and a slowly-increased concentration depen- 
dency take place up to fields E ~20 kV/cm in agreement 
with the experimental data. [7j This is because an in- 
terplay between the emission and absorption of optical 
phonons (the last process is proportional to the nonzero 
Planck number of T- and if- modes). For lower tem- 
peratures, I(V) tends to be saturated (starting E ~20 
V/cm at T <150 K) and a negative differential conduc- 
tivity (NDC) appears in the transition region {E ~ 5 - 
15 V/cm at 77 K). Since a predominant recombination 
of high-energy carriers due to optical phonon emission 
at low temperatures, a depletion of concentration takes 
place below 250 K. The semi-insulating regime of conduc- 
tivity with a residual concentration about 10 s cm~ 2 , an 
energy per particle less 80 meV, and with a low saturated 
current is realized at 77 K and E ^20 V/cm. Beside of 
the peculiarities of the steady-state response, a two-scale 
temporal evolution takes place due to a fast relaxation 
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via optical phonons and a slow relaxation via acoustic 
phononsand thermal radiation. A scale of transient pro- 
cess is strongly dependent on temperature and on fields 
applied. At low temperatures a steady-state distribution 
is established during ~10 ns time scale and the transient 
response can be measured directly. 

The paper is organized as follows. The basic equations 
governing the heating of carriers under an abrupt switch- 
ing off a dc electric field are considered in Sec. II. Tempo- 
ral evolution of nonequilibrium distributions for different 
temperatures and electric fields is described in Sec. III. 
Description of the temporal and steady-state responses, 
including the nonequilibrium concentration and energy as 
well as the current-voltage characteristics, is presented in 
Sec. IV. The concluding remarks and discussion of the 
assumptions used are given in the last section. 



II. KINETIC APPROACH 

Nonequilibrium electrons and holes in an intrinsic 
graphene are described by coinciding distributions f p t 
because their energy spectra are symmetric and the scat- 
tering mechanisms are identical. [9( An evolution of these 
distributions under a homogeneous electric field E t is 
governed by the quasiclassical kinetic equation 
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Here the collision integrals Jj(ft\p) describe the relax- 
ation of carriers caused by the elastic scattering on struc- 
ture disorder (j = d), the intraband quasielastic scatter- 
ing on acoustic phonons (j — qe), and the intra- and 
interband emission and absorption of optical phonons 
(j = intra, inter). In addition, the contribution J r (ft\p) 
describes the generation-recombination processes due to 
interband transitions induced by thermal radiation. At 
t < 0, we use the equilibrium initial condition f p t<o = 

f p eq ^ at temperature T where f p = {ex.p(p/px) + l} -1 
and pt — T/v is written through the characteristic ve- 
locity v ~ 10 8 cm~ 2 . 

For typical graphene structures, the momentum re- 
laxation dominants over the other processes listed, i.e. 
Vd 3> Vj, j 7^ d where Vj stands for the relaxation 
rate of the j-th scattering channel. Under the condition 
eEr m <C p, where r m ~ v^ 1 is the momentum relaxation 
time and vp is the energy of hot carriers, the anisotropy 
of distribution is weak, f pt ~ f pt + A/ pt . Here we sepa- 
rated the isotropic part of distribution, f pt , and a weak 
anisotropic contribution, A/ pt . Within the local time ap- 
proximation, when E t increase slowly over r m -scale, we 
obtain the asymmetric part of distribution at t > as 



The phenomenological rate v p is written here through the 
characteristic velocities Vd and vq for the case of the com- 
bined elastic scattering by short- and long-range disorder, 
see Ref. 12 for details. For the model of the Gaussian 
disorder with the correlation length l c we use the dimcn- 
sionless function *&(z) — e~ z I\{z 2 )jz 2 with the first- 
order Bessel function of an imaginary argument, I\(z). 

Performing the averaging over in-plane angle (such av- 
eraging symbolized by overline) in kinetic equation (1) 
and neglecting the weak contribution of A/ pt in the right- 
hand side, one obtains the kinetic equation for symmetric 
distribution f pt as follows 
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Because the elastic scattering does not affect the sym- 
metric distribution due to the energy conservation law, 
the only non-elastic mechanisms (j ^ d) are responsi- 
ble for relaxation of f p t- The Joule heating contribution 
(oc E 2 term in the left-hand side) is expressed through 
the asymmetric correction (2) and it can be transformed 
into 
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Eqs. (3, 4) should be solved at t > with the initial con- 
dition / p t=o = fp eq ^ and the use of the collision integrals 
described below. 

Within the quasielastic approximation, [H, [l3| the en- 
ergy relaxation via acoustics phonons is described by the 
Fokker-Planck collision integral given by 
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Here the rate v p = (s/v) 2 v ac p/Ti is written through 
the sound velocity s (v/s ~137) and the characteristic 
velocity v ac ~2.4xl0 5 cm/s at the nitrogen temperature, 
moreover v ac oc T. The interband transitions caused by 
the thermal radiation are described by the collision inte- 
gral 

(2 ■ 



Mft\p) = v v U) [N 2p/pT [l - 2f pt ) - f 2 t ] 



(6) 



where the rate of spontaneous radiative transitions, 
Vp = v r p/%, is written through the characteristic ve- 
locity v r ~41.6 cm/s for the case of graphene on the 
SiC>2 substrate. 

The scattering by optical phonons is described by the 
collision integral written through the intra- and inter- 
band parts, Jintra(ft\p) + J inter (ft\p), as follows (see eval- 
uation in Refs. 14 and 15) 

Jintra{ft\p) = ^ [(N n + 1)^+^(1 ~ fpt)fp+p n t 



eE ■ p / df p t \ v =??^Pq,(Ph\ 
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and 
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Here N v = [exp(huj v /T) — is the Planck distribution 
of 77th phonon mode with the energy fkjj v at tempera- 
ture T and p v = fiu^/v is the characteristic momentum. 
Summation over 77 involves both the intra- and interval- 
Icy transitions (marked by r\ — T and K respectively) 
taking into account the zone-center and zone-boundary 
phonon modes. Since the averaging over in-plane an- 
gle in the general form of the collision integral, [l4], EH] 
one obtains the relaxation rates up 1 ' ~ 9(p)v ri p/h and 

Vp ~ 9(p)v n p/h, which are proportional to the den- 
sity of states. The ^-function here allows the interband 
transitions, if p < p v (in the passive region). Below we 
suppose the same relaxation rates for the intra- and inter- 
valley transitions, v v ~ v n and we use the characteristic 
velocities «r — 10 6 cm/s and vk — 2 x 10 6 cm/s which 
are in agreement of the previous results. [H], EH 

The concentration of electrons (holes) , n t , and the en- 
ergy and current densities, £t and It, are determined 
through f pt and A/ p t according to the standard rela- 
tions: 
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The factor 4 here takes into account the spin and valley 
degeneracy and v p = vp/p is the carrier velocity. Using 
the local time approach, see Eq. (2), and introducing 
the nonlinear conductivity at according to It = crtE t , 
one obtains 
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For the case of the short-range scattering (l c — >0 and Vd = 
v c i + vq), when v v ~ vdp/fi, this integral is transformed 
into a t ~ (e 2 v/Trh 2 Vd)f P =ot, i.e. a t oc ft for the low- 
energy region [vp <60 - 100 meV, see Fig. 1(a) below]. 



III. TRANSIENT EVOLUTION OF 
DISTRIBUTION 



First, we consider numerical solution of the Cauchy 
problem for the nonlinear differential equation with the 
finite-difference terms given by Eqs. (4)- (8). Using the 
iteration scheme [l6| at t > and the equilibrium initial 
condition one obtains the transient distribution f pt for 
different electric fields and temperatures. Under these 
calculations we use the momentum relaxation rate given 
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FIG. 1: (Color online) (a) Momentum relaxation rate v v ver- 
sus energy vp for l c —5 nm (1) 7.5 nm (2), and 12.5 nm (3). 
(b) Conductivity versus gate voltage for the same l c as in 
panel (a). Squares are experimental points from Fig. 2a of 
Ref. 7. 



by Eq. (2) at different correlation lengths, l c =5 - 12.5 
nm, with the characteristic velocity Vd — 2.6 x 10 7 cm/s, 
which is correspondent to the maximal sheet resistance 
~3.3 kO and Vo/vd ~0.035. In Figs. 1(a) and 1(b) we 
plot the relaxation rate v v versus energy vp and the linear 
conductivity a versus gate voltage, V g , for l c —5, 7.5 
and 12.5 nm. As it is shown in Fig. 1(b), the curve 
2 is in agreement with the experimental data of Ref. 7 
and below we use these parameters for description of the 
nonlinear response. 

Transient evolution of distribution f pt , from equilib- 

p( e <?) 



(9) rium form fp to steady-state function f p = f } 



pt— ^001 is 

shown in Figs. 2(a) - 2(f) for different temperatures and 
electric fields. At room temperature, the transient pro- 
cess takes place during 1 - 10 ps time scale and this pro- 
cess becomes slower for higher fields because an interplay 
between emission and absorption of optical phonons. At 
lower temperatures, the evolution times are longer, about 
1 ns at 150 K and 50 ns at 77 K, but these time scales 
become shorter for higher fields in contrast to the room 
temperature case. It is because N v < 10~ 5 at T <200 K 
and the absorption of optical phonons is negligible. The 
maximal distribution f p= ot decreases during the transient 
process, moreover this effect is enhanced for higher fields 
and lower temperatures. 

The steady-state distribution functions, f p , which are 
established after the transient process, are shown in Fig. 
3 for different field strengths and l c =7.5 nm at tempera- 
tures 300 K, 150 K, and 77 K. In contrast to the exponen- 
tial decay of the equilibrium distribution at E = shown 
as the dotted lines, at E ^ one obtains a smeared dis- 
tribution with a low maximal value / p =o < 10~ 2 . At 
room temperature, carriers are distributed over energies 
3> ftZJri, if E > 2 kV/cm [see Fig. 3(a)] due to competi- 
tion between absorption and emission of optical phonons. 
Here and below, hZJ^ stands for the lowest optical phonon 
energy. At low temperatures, the distributions are lo- 
cated in the passive region of energies vp < ?HU V [see 
Figs. 3(b) and 3(c)] due to effective emission of optical 
phonons. 
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FIG. 2: (Color online) Distribution function f p t versus energy 
vp for times marked at different temperatures and electric 
fields: (a) T =300 K and E =6 kV/cm; (b) T =300 K and 
E =12 kV/cm; (c) T =150 K and E =60 V/cm; (d) T =150 K 
and E =120 V/cm; (e) T =77 K and E =10 V/cm; (f) T =77 
K and and E =20 V/cm. Dotted and solid black curves are 
correspondent to initial and final distribution, fp Eq ^ and f p , 
respectively. 
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FIG. 3: (Color online) (a) Distribution function f p versus vp 
for T =300 K at different electric fields (marked), (b) The 
same for T =150 K. (c) The same for T =77 K. Dotted lines 
are correspondent to the equilibrium distribution. 
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FIG. 4: (Color online) Transient evolution of concentrations 
nt (normalized to their equilibrium values, nr) for l c =7.5 nm 
and different electric fields (marked) at T =300 K (a), 150 K 
(b), and 77 K (c). 



IV. RESULTS 

The nonlinear responses introduced by Eqs. (9, 10) are 
analyzed below. Both the transient characteristics and 
the steady-state dependencies of concentration, energy 
and current densities on temperature and on dc field are 
described. 



A. Concentration and energy vs t and E 

Transient evolution of concentrations n t , normalized 
to the equilibrium concentration in an intrinsic graphene 



[n T ~ 0.52(T/hv) 2 varies between 8.1 x 10 10 cm" 2 and 
5.4 x 10 9 cm~ 2 for T =300 - 77 K] is shown in Figs. 4(a) 
- 4(c) for the same temperatures and electric fields as in 
Fig. 2. At room temperature, n t increases with t and the 
saturation regime is realized at longer times for higher 
fields (t >0.1 ns for E =12 kV/cm). It is because the 
distribution is moved away from the passive region [see 
Figs. (2a) and (2b)] where the generation-recombination 
processes take place. At lower temperatures, n t decreases 
during temporal evolution and the saturation regime is 
realized for time intervals >0.5 ns or >50 ns at T =150 K 
or 77 K, respectively. The depletion of concentration is 
realized because an effective recombination via the opti- 
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cal phonon emission appears at vp > huJ v /2. All tempo- 
ral dependencies appear to be faster during initial stages 
of evolution at high fields, i. e. a two-stage transient pro- 
cess takes place because an increasing of the relaxation 
rates with energy. 
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FIG. 5: (Color online) (a) Concentrations of carriers n (nor- 
malized to their equilibrium values, nr) versus electric field 
E for T =300 K (a), 225 K (b), 150 K (c), and 77 K (d). 
Curves 1-3 are correspondent to l c =5, 7.5, and 12.5 nm. 

Variation of the normalized nonequilibrium concentra- 
tion n/riT with field applied is shown in Figs. 5(a)-(d) 
for different temperatures and l c . At room temperature 
n increases with E due to the thermogeneration process 
(first term of Eq. (8) with N K ~ f .4 x 10 -3 ). At lower 
temperatures, if T <250 K when N v < 1CT 4 , the ther- 
mogeneration is negligible and Jinter < 0. As a result, 
depletion of concentration takes place as E increases. In 
a high-field region [at E >50 V/cm if T =77 K, see Fig. 
5(d)], one arrives to the saturation regime with concen- 
trations lower 10 8 cm~ 2 . The generation-recombination 
processes appears to be more effective if l c >10 nm, so 
that n decreases faster if E increases. 

It is convenient to characterize the energy distribution 
by the averaged energy per carrier, £t/nt, normalized to 
their equilibrium value £t/tit, i- e. we consider below 
the ratio 
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(11) 



where Et/tlt — 0.22T. Transient evolution of £t is plot- 
ted in Figs. 6(a) - 6(c) for the same T and E as in Figs. 2 
and 4. Notice, that the energy per carrier saturates faster 
in comparison to concentration, c. f. Figs. 6 and 4. The 
steady-state ratio £ = £t->oo is shown in Fig. 7 versus 




FIG. 6: (Color online) Transient evolution of energy per car- 
rier normalized to their equilibrium values, £t given by Eq. 
(11) for l c =7.5 nm and different electric fields (marked) at 
T =300 K (a), 150 K (b), and 77 K (c). 
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FIG. 7: (Color online) The ratio £ = £ t ^oo given by Eq. (11) 
versus electric field E for T =300 K (a), 150 K (b) and 77 K 
(c). Curves 1-3 are correspondent to l c =5, 7.5, and 12.5 
nm. 



electric field for different l c and temperatures. After a 
fast increasing of £ with E at low fields, one obtains a 
slow dependency of the energy per carrier versus E. The 
dependencies on l c are similar to ones shown in Fig. 5. 



B. Current-voltage characteristics 

We turn now to consideration of the current-voltage 
characteristics, I t versus E t normalized to length of sam- 
ple. Transient evolution of current I t is shown in Fig. 8 
for the same T and E as in Figs. 2, 4, and 6. Once again, 
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FIG. 8: (Color online) Transient evolution of current- volt age 
characteristics for the same conditions as in Fig. 4. 
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FIG. 9: (Color online) Current-voltage characteristics for the 
same conditions as in Fig. 5. Experimental points from Fig. 
2(b) of Ref. 7 are shown in panel (a) as squares. 



ps, c. f. Figs. 4(a) and 8(a). 

The steady-state current- voltage characteristics, 1(E), 
are shown in Fig. 9 under the same conditions as in Figs. 
5 and 7. For the low-field region (about tens V/cm at 
T >200 K and -V/cm at lower T), the 1(E) dependen- 
cies are in agreement with the results reported in Ref. 9. 
In the high-field region, the second Ohmic law with the 
resistance ~3 kJ7 at T =300 K takes place, while the sat- 
uration of current with I <10 /zA/cm takes place at low 
temperatures (a semi-insulating regime of transport). In 
an intermediate field region (—10 - 30 V/cm at 150 K or 
—5 - 15 V/cm at 77 K), a decreasing 1(E) characteristics 
are realized, i.e. the NDC regime of response takes place. 
Both the effective conductivity in the second Ohmic law 
region and the differential conductivity under the NDC 
regime are increased with l c . 
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FIG. 10: (Color online) Differential conductivities dl /dE at 
temperatures T =150 and 77 K (dashed and solid curves) 
for l c =5, 7.5, and 12.5 nm (curves 1-3, respectively). Inset 
shows region of parameters, E and T (shaded) for which NDC 
regime takes place at l c =7.5 nm. 

In Fig. 10 we plot the differential conductivity dl /dE 
versus E for the temperatures 150 and 77 K when the 
NDC regime takes place. Notice that dl/dE shows a 
strong dependency on l c in the low-temperature region 
[c. f. Figs. 9(c) and 9(d)]. The NDC region over (E, 
T)-plane, where (dl/dE) < is shaded in the inset for 
l c =7.5 nm. A possibility for the development of a spatial 
instability under NDC regime requires a special investi- 
gation similar to the bulk case. [17| But due to a small 
values of \dI/dE\ (<100 fiS at T >20 K) one may expect 
that an inhomogeneous distribution is not developed for 
the typical size of samples. 



one obtains the two-stage evolution over time scales sim- 
ilar to those for the n t and £ t dependencies. At low tem- 
peratures (150 K and 77 K) I t decreases monotonically, 
similarly to n t . At room temperature, It decreases fast 
during an initial stage of the transient process because en- 
ergy and momentum relaxation rates increase. For longer 
times, It increases similarly to the nt dependency at t >5 



V. CONCLUSIONS 

Summarizing the consideration performed, we have ex- 
amined the heating of high-energy carriers taking into 
account the effective intra- and interband transitions 
caused by the optical phonons. The interplay of Joule 
heating and recombination of electron-hole pairs gives 
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rice to (i) the depletion of carriers' concentration and 
to (ii) the negative differential conductivity at low tem- 
peratures. The model developed explains the nearly lin- 
ear current-voltage characteristic at room temperature, 
in agreement with the experimental data for intrinsic 
grapnene. Q Beside of this, temporal evolution of re- 
sponse under an abrupt switching off a dc field lasts over 
a nanosecond time scale, so that a direct measurement of 
transient response open up the way to verify the relax- 
ation and recombination mechanisms. 

Further, we compare the consideration presented with 
the other results published in order to stress that the non- 
linear response is sensitive to a relative contribution of 
different relaxation processes and to a geometry of mea- 
surements (size of sample, contacts) . The last regime was 
considered in short structures, [3|-|5[ see a recent data in 
Ref. 18 and a description of the ballistic limit of trans- 
port was performed in Refs. 19. In the case of monopolar 
transport @ the concentration is fixed and I(V) charac- 
teristics are only determined by the momentum and en- 
ergy relaxation mechanisms, as it was discussed in Refs. 
10 and 11. Last but not least, the Joule heating is de- 
termined by the momentum scattering mechanism and it 
should be verified from the linear conductivity measure- 
ments, see Fig. 1 and Refs. 7 and 12. An additional 
uncertainty appears due to the use of any model of mo- 
mentum relaxation because a microscopic mechanism re- 
mains under debates. 0, Ho| 

Next, we discuss the assumptions used. The local time 
approach and the weak anisotropy of distribution used 
in (2) are valid under the dominant momentum scat- 



tering. We omitted the scattering by surface phonons 
of the substrate in agreement with the consideration of 
the experimental data. HJ If such a contributions 
essential, it can shift the peculiarities considered (deple- 
tion of concentration and negative differential conductiv- 
ity) to lower temperatures because the optical phonon 
energy ^55 meV for the SiC>2 substrate. The carrier- 
carrier scattering is unessential here because the distri- 
butions obtained are spreaded over a range of energies up 
to hur,K, so that f p <C 1. Also one can neglect by the 
long-range disorder [22] because there is no low-energy 
particles if / p =o *C 1. Since Vd 3> vr,K> one can neglect 
the scattering by optical phonons in comparison by the 
elastic scattering under consideration of the momentum 
scattering. In addition, a heating of phonons is neglected 
under the condition of effective phonon thermalization. 
['2'.'l Besides, we are restricted by the quasiclassical ap- 
proach, neglecting a mix of bands due to a strong field 
(see [18| and the quasiclassical condition in Appendix of 
Ref. 9a), which can be essential in a lateral pn-junction. 
The simplifications listed do not change either the char- 
acter of the high-field distributions or the numerical es- 
timates for the current-voltage characteristics. 

To conclude, both a detail theoretical consideration, 
including numerical modeling for different scattering and 
recombination mechanisms, and an experimental study 
of heating at low temperatures are timely now. An ob- 
servation of the peculiarities obtained and a verification 
of relaxation processes can be performed for a homoge- 
neous samples with the four-point geometry of contacts. 
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